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Abstract 

The Analytic Hierarchy Process (AHP) is widely used for decision making in- 
volving multiple criteria. Eisner and van den Driessche |1Q|, [TT] introduced a max- 
algebraic approach to the single criterion AHP. We extend this to the multi-criteria 
AHP, by considering multi-objective generalisations of the single objective optimi- 
sation problem solved in these earlier papers. We relate the existence of globally 
optimal solutions to the commutativity properties of the associated matrices; we 
relate min-max optimal solutions to the generalised spectral radius; and we prove 
that Pareto optimal solutions are guaranteed to exist. 
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1 Introduction 



The analytic hierarchy process (AHP) is a method for ranking alternatives in multi- 
criteria decision making problems. Developed by Saaty [26J, it consists of a three layer 
hierarchical structure: the overall goal is at the top; the criteria are in the next level; and 
the alternatives are in the bottom level. The AHP has been used in many different areas 



"Hamilton Institute, National University of Ireland, Maynooth, Co. Kildare, Ireland. Email: 
buket. benek@nuim.ie 

t Corresponding author. Hamilton Institute, National University of Ireland, Maynooth, Co. Kildare, 
Ireland. Email: oliver.mason@nuim.ie 

■^University of Birmingham, School of Mathematics, Edgbaston B15 2TT. Email: sergiej@gmail.com 



1 



including manufacturing systems, finance, politics, education, business and industry; for 
more details on the method, see the monographs by Saaty- Vargas and Vaidya-Kumar J27J 



The essence of the AHP can be described as follows. Given n alternatives we construct 
a pairwise comparison matrix (PC-matrix), A > for each criterion, in which indicates 
the strength of alternative % relative to alternative j for that criterion. A PC-matrix with 
the property that a^aji = 1 for all i ^ j and an = 1 for all i is called a symmetrically 
reciprocal matrix (SP-matrix) [12J. (Note that this abbreviation might clash with the 
strongly regular matrices of Butkovic [5], but not in this paper.) 

Once an SP-matrix A is constructed, the next step in the AHP is to derive a vector 
(wi, . . . , w n ) of positive weights, which can be used to rank the alternatives, with Wi 
quantifying the weight of alternative %. As observed by Eisner and van den Driessche [TO] , 
the ideal situation is where a^- = Wi/wj, in which case the S'P-matrix is transitive. In 
practice, this will rarely be the case and it is necessary to approximate A with a transitive 
matrix T, where Uj = Wi/wj for some positive weight vector w = (wi, . . . ,w n ). The 
problem is then how to construct T given A. Several approaches have been proposed 
including Saaty's suggestion to take w to be the Perron vector of A, or the approach of 
Farkas et al. [12], which chooses w to minimise the Euclidean error Y2{ a ij~ w i/ W j) 2 - Eisner 



and van den Driessche [TDl ITT] suggested selecting w to be the max algebraic eigenvector of 
A. This is similar in spirit to Saaty's approach and also generates a transitive matrix that 
minimises the maximal relative error max |a^- — Wi/wAj ' a^. As noted in [TTJ, minimising 
this functional is equivalent to minimising 



The different approaches to approximating an SP-matrix A with a transitive matrix T 
will in general produce different rankings of the alternatives. The question of how these 
rankings are affected by the choice of scheme is considered in the recent paper of Ngoc [29] . 

In the classical AHP involving multiple criteria, a set of SP-matrices is constructed: 
one for each criterion. One additional S'P-matrix is constructed based on comparisons 
of the different criteria. Once weight vectors are obtained for each individual criterion, 
these are then combined using the entries of the weight vector for the criteria-comparison 
matrix. As an illustration, we take the following numerical example from Saaty [26] and 
show how the Perron vectors of the comparison matrices are used to construct a weight 
vector. 



30J. 




l<i,j<n 



(i) 
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Example 1.1. The problem considered is deciding where to go for a one week vacation 
among the alternatives: 1. Short trips, 2. Quebec, 3. Denver, 4- California. Five criteria 
are considered: 1. cost of the trip, 2. sight-seeing opportunities, 3. entertainment, 4- 
means of travel and 5. dining. The PC-matrix for the criteria and its Perron vector are 
given by 

" 1 1/5 1/5 1 1/3 1 I" 0.179 

5 1 1/5 1/5 1 0.239 
C= 5 5 11/5 1 and c= 0.431 

1 5 5 1 5 0.818 
3 1 11/5 1 0.237 

The above matrix C describes the pairwise comparisons between the different criteria. For 
instance, as C21 = 5 ; criterion 2 is rated more important than criterion 1; C32 = 5 indicates 
that criterion 3 is rated more important than criterion 2 and so on. The vector c contains 
the weights of the criteria; in this method, criterion 4 is given most weight, followed by 
criterion 3 and so on. 

The SR-matrices, A ± , A 5 , for each of the 5 criteria, their Perron vectors and cor- 
responding ranking schemes are given below. For instance, for criterion 1, the first alter- 
native is preferred to the second as the (1,2) entry of A\ is 3. Similarly, for criterion 3, 
the 4th alternative is preferred to the 1st as the (4, 1) entry of A 3 is 2. 

For the cost of the trip: 



A, 



1 3 
1/3 1 
1/7 1/6 



7 
6 
1 



9 
7 
3 



_ 1/9 1/7 1/3 1 

For the sight-seeing opportunities: 



A, 



1 1/5 1/6 1/4 

5 12 4 

6 1/2 1 6 
_4 1/4 1/6 1 





' 0.877 " 




0.46 




0.123 




. 0.064 _ 




" 0.091 " 




0.748 




0.628 




0.196 



1 > 2 > 3 > 4 



2 > 3 > 4 > 1 



For the entertainment: 



A, 



1 7 7 1/2 
1/7 1 1 1/7 
1/7 1 1 1/7 

2 7 7 1 





0.57 " 




0.096 




0.096 




0.81 



4 > 1 > 2 = 3 
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For the means of travel: 





1 


4 


1/4 


1/3" 




' 0.396 _ 


A A = 


1/4 


1 


1/2 


3 




0.355 




4 


2 


1 


3 




0.768 




3 


1/3 


1/3 


1 




. 0.357 . 



For the dining: 





1 


1 


7 


4 




' 0.723 _ 




1 


1 


6 


3 




0.642 




1/7 


1/6 


1 


1/4 




0.088 




. 1/4 


1/3 


4 


1 




_ 0.242 _ 



3 > 1 > 4 > 2 



1 > 2 > 4 > 3 



To obtain the overall weight vector, we compute the weighted sum ^2 (kv^ 1 ' . This gives 

i=l 



w 



0.919 
0.745 
0.862 
0.757 



with the associated ranking: 1 > 3 > 4 > 2. 



Our work here is inspired by the max-algebraic approach to the AHP introduced by 
Eisner and van den Driessche pJJ, [UJ and extends it in the following manner. In [TTJ1 [TTj . 
the max eigenvector is used as a weight vector for a single criterion and it is shown 
to be optimal in the sense of minimising the maximal relative error as discussed above. 
This work naturally raises the question of how to treat multiple criteria within the max- 
algebraic framework. We address this question here by considering the multi-criteria AHP 
as a multi-objective optimisation problem, in which we have an objective function of the 
form (pQ) for each criterion (and associated 572-matrix) . Rather than combining individual 
weight vectors as in Example 11.11 we consider three approaches within the framework of 
multi-objective optimisation, and use the optimal solution as a weight vector in each case. 
The advantage of this approach is that the weight vector can be interpreted in terms of the 
maximal relative error functions (CEJ) associated with the 57?-matrices given as data for the 
problem. The optimisation problems we consider are the following. First, we investigate 
the existence of a single transitive matrix with a minimum distance to all matrices in the 
set simultaneously. We remark that this amounts to finding a common subeigenvector of 
the given matrices. Clearly, this will not in general be possible. The second problem we 
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consider is to obtain a transitive matrix that minimises the maximal distance to any of 
the given ^-matrices. The third problem concerns the existence of a transitive matrix 
that is Pareto optimal for the given set of matrices. To illustrate our results, we revisit 
Example 11.11 towards the end of the paper. 

2 Notation and Mathematical Background 

The set of all nonnegative real numbers is denoted by M + ; the set of all n-tuples of non- 
negative real numbers is denoted by R™ and the set of all n x n matrices with nonnegative 
real entries is denoted by IR" xn . We denote the set of all n-tuples of positive real numbers 
by int(lR™). For A e M+ Xri and 1 < i,j < n, a {j refers to the (i, j) th entry of A. The 
matrix A = [ay] is nonnegative (positive) if > (a^ > 0) for 1 < i,j ' < n. This is 
denoted by A e M+ Xn (A > 0). 

The weighted directed graph of A is denoted by D(A). It is an ordered pair (N(A), E(A)) 
where N(A) is a finite set of nodes {1,2, ...,n} and E(A) is a set of directed edges, with 
an edge (i, j) from i to j if and only if a»j > 0. A path is a sequence of distinct nodes 
ix, i 2 , ■ ■ ■ , ik of length k — 1 with the weight a ili2 a i2 i 3 ■ ■ ■ a ik _ lik , where (i p , i p+ \) is an edge 
in D(A) for p = 1, . . . , k — 1. It is standard that A is an irreducible matrix if and only if 
there is a directed path between any two nodes in D(A). A cycle T of length k is a closed 
path of the form ii, «2, iki H- The k th root of its weight is called its cycle geometric 
mean. For a matrix A 6 R™ x ™, the maximal cycle geometric mean over all possible cycles 
in D(A) is denoted by fi{A). A cycle with the maximum cycle geometric mean is called 
a critical cycle. Nodes that lie on some critical cycle are known as critical nodes and de- 
noted by N C '(A). The set of edges belonging to critical cycles are said to be critical edges 
and denoted by E C (A). The critical matrix of A [HIE], A c , is formed from the submatrix 
of A consisting of the rows and columns corresponding to critical nodes as follows. Set 
a ?j = a ij ^ ihj) h es on a critical cycle and af^ = otherwise. We use the notation D C (A) 
for the critical graph where D C (A) = D(A C ) = (N C (A),E C (A)). 

The max algebra consists of the set of nonnegative numbers together with the two basic 
operations a © b = max(a, b) and a £g> 6 = ab. This is isomorphic to the max-plus algebra 
[U, [5] via the natural isomorphism x — > log(x). These operations extend to nonnegative 
matrices and vectors in the obvious manner pQ [2j [7j. For A in IR" xn , the eigenequation 
in the max algebra is given by A <g> x = Xx, x > 0, A > 0. fi(A) is the largest max 
eigenvalue of A [5] . If A is irreducible, then it is the unique max eigenvalue of A and there 
is a positive max eigenvector x G int(M") corresponding to it [H [2]. The eigenvector v is 
unique up to a scalar multiple if and only if A c is irreducible. 
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Observe that an 572-matrix A is irreducible and fi(A) > 1 [10J. Although our primary 
interest is in 572-matrices, it is noteworthy that many of our results also hold true for 
non-zero reducible matrices. 

For A E M+ Xn with /jl(A) < 1, I ® A® A® ® ... converges to a finite matrix called the 
Kleene star of A given by A* = I © A © A% © ... © A 11 ' 1 where fi(A*) = 1 [D [28] . Here, 
a*^ is the maximum weight of a path from % to j of any length [16j (if i ^ j), and 
denotes the k th max-algebraic power of A. Note that for each A E IR" xn , if A* is finite 
then the max-algebraic sum of all of the columns of A* is positive. For A E IR™ xn , the set 
of subeigenvectors of A associated with fi(A) is called a subeigencone of A and denoted 
by V*(A) = {y e \ A (g) y < ^{A)y} [28J . It was shown in Proposition 2.5 of [28] 
that for A E M r | xn with fj,(A) = 1, V*(A) = V(A*) = span e (A*) where V(A*) denotes the 
eigencone of A* consisting of its max eigenvectors. Note that the above-mentioned max- 
algebraic sum of all columns of A* is in V*(A) [28], so V*(A) contains positive vectors. 
Note that if /i(A) > 0, we can normalise A by fi(A) and V*(-^) = V*(A). 

To the authors' knowledge, max-algebraic subeigenvectors appeared in the works of 
Gaubert [T3l H3] . However, they can be traced back to earlier works on nonnegative 
matrix scaling, see references in Butkovic-Schneider [6]. 

For A £ ]R™ xn we will consider the following set, which was introduced in [11] 



For the special case of r = fi(A), Ca^(a) is denoted by Ca [II]- Obviously Ca is the 
positive part of V*(A) (which is non-empty as we argued above), and it coincides with 
V*(A) when A is irreducible. To be consistent with the notation of [TT] . we recall the 
definition of the normalised set 



As above, T>a is used to denote the special case where r = /jl(A). 

The relations between the sets Ca,t-, the error function and n{A) were clarified by 
Eisner and van den Driessche [H] and are recalled in the following propositions, which 
we easily extend (based on [5j El [28]) to the general reducible case. 

Proposition 2.1 (cf. [H] Lemma 2). Let A £ M" xri be nonzero. Then: 



C A>r = {x £ int(R+) | A © x < rx}. 



(2) 



T> A ,r = {X £ C A ,r | X\ = !}• 



(3) 



(i) C A , r ^ 



r > 0, r > /i(A); 
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Proof. To prove(i), exploit |5] Theorem 1.6.29 stating that 

/j(A) = min{A \ A® x < \x, x G int(R" ), }, (4) 

when ii{A) > (based on Butkovic-Schneider [6j Theorem 2.6). In the trivial case r = 
n{A) = 0, we have V*(A) = {x \ A £g> x = 0}, which consists of all vectors x such that 
X{ ^ if and only if the 2 th column of A is zero. In this case Ca,t = T^a,t = 0, unless 
A = (which we exclude). 

The result of (ii) follows from the definitions of ca{x) and C^, r given in ([1]) and ([2]). □ 

See also Gaubert [13] Ch. IV Lemma 1.3.8, Krivulin [2D] Lemma 1, and Nussbaum [21] 
Theorem 3.1 (in a more general nonlinear context) for closely related statements. 

Proposition 2.2 (cf. [11], Theorem 1 part 6). Let A G M™ x ". Then x G V>a is unique if 
and only if A c is irreducible and N C '(A) = N(A). 

Proof. Proposition 12.11 (i) implies that we have T>a = Ca = when n{A) = 0, so we can 
assume (i(A) > and, further, fi(A) = 1. According to [2H] Theorem 2.8, 

V*(A) = V{A*) = \ Atf*© X jg j lA^GR+l. (5) 

i£M(A) i£N c (A) ) 

Here g % is the i th column of A*. The subset M(A) C {1, . . . ,n} is such that for each 
(maximal) irreducible submatrix of A° there is a unique index of that submatrix in 
M(A). By the same theorem of [25] , based on the classical results in [TJ [7j, columns of 
A* with indices in the same irreducible submatrix of A c are proportional to each other, 
and there is no proportionality between different g l appearing in ((5]) (moreover, these g l 
are strongly linearly independent [5]). 

'If" : As A c is irreducible and all nodes are critical, all columns of A* are proportional 
to each other, and ([5]) shows that V*(A) is just one ray. A > A c is irreducible as well, so 
V*(A) = Ca and V>a is the unique vector on that ray with x\ = 1. 

n 

"Only if": Let A c be reducible, or let N C (A) ^ N(A). The vector z = g i is 

i=i 

positive. Consider g % © ez for all i, taking small enough e so that all entries of ez are 
less than any nonzero entry of A*. In this case the positive entries of g % do not change, 
and if g l and gi are not proportional then neither are g % © ez and gi © ez. After suitably 
normalising these vectors, we obtain two different vectors in T>a- □ 

Next, let \1/ C M" xn be a finite set of nonnegative matrices given by 

^ = {A 1 ,A 2 ,...,A m }, 3i:Ai^Q (6) 
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Given let ^ denote the set of all products of matrices from \I/ of length p > 1. 
Formally, ^ = {A^ ® • • • © A, p I < jk < m for 1 < k < p}. Using this, the max version 
of the generalised spectral radius [2TJ [25] is defined by 

/t(^) = limsup(max fi(ip))p. (7) 

Before stating the next theorem, let S be the matrix given by 

S = 0A = A x © A 2 © ... © A m . (8) 

Note that S 1 > if at least one Aj > 0. Moreover, fi(S) > if at least one /i(Aj) > 
and S is irreducible if at least one Ai is irreducible: these represent the main cases in 
which we are interested. 

Theorem 2.1. Let ^ be given by (0|) and S be given by Then, /t(\I / ) = n(S) . ( Gaubert fL^ , 
Benek Gursoy and Mason JEjl) 

Inspired by the approach to the single criterion AHP adopted in [11], we associate a 
set &!J tr with the set \I/ of nonnegative matrices and show how the geometric properties 
of Ca,t discussed in [11] extend to this new setting. 

Define 

m 

C*, r = {xe int(R+) I e Ai (x) < r for all A; G ^} = f] C Ai>r . (9) 

i=i 

We also consider the set of normalised vectors: 

£>*, r = {x G C*, r | X! = 1}. (10) 

We will use the notations Cy for C^^) and 2^ for X^,^). The following result shows 
that the set C^ jr and the set Cs, r ar e equal. This will allow us to readily extend properties 
of Ca,t established in [11] to sets of matrices. 

Theorem 2.2. Consider the set \1/ C IR" xn given by and let S be given by §B\). Then: 

C^ t r = Cs,r- 

Proof, (i): Let x in C^ <r be given. Then, e^(x) < r which implies Aj © x < rx for each 
Aj G ^ from (ii) in Proposition 12.11 Taking the maximum of both sides from 1 to m, we 
see that S © x < rx. It follows that x G C$ jr - Thus, C\& )r C Cs, r - 

Now choose some x G Cs, r - Then es(a;) < r from (ii) in Proposition 12.11 Since 
£Ai(x) < es(%) < r for all 1 < i < m, we obtain x G C^. r . Thus, Cs, r C C$, r . Hence 
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The following corollary extends Proposition 12.11 to a set of nonnegative matrices 
Since frfif) = n{S) by Theorem 12.11 (ii), the corollary is an immediate consequence of 
Theorem 12.21 

Corollary 2.1. Consider the set ^ C M" xn given by fifij), and let S be given by (dp. Then: 

(i) C*, r ^ r > 0, r > 

(ii) x G C^ >r •<=>- x G int(M"), eg (a;) < r. 

Theorem 12.21 establishes that C%/ jT = Cs, r - It is immediate that we also have P* jr = 
Ds,n Cq, = Cs and Vq, = T>s- Therefore, studying these sets for a collection of nonnegative 
matrices reduces to studying the sets associated with the single matrix S. This fact means 
that the properties of Ca,t discussed in Theorem 1 of [11] can be directly extended to C$ jr . 
We state some of these in the following theorem, noting that property (i) does not require 
irreducibility. 

Theorem 2.3 (cf. [llj Theorem 1). Let C R™ x " be given in (0|), and let S be given by 
(TJj. Then: 

(i) &i, tr and T>^ r are convex and max-convex; 

(ii) consists of only one vector if and only if S c is irreducible and N C (S) = N(S). 
(Hi) If all matrices in $ are irreducible then P* >r is compact. 

Proof, (i): (see [28] Proposition 3.1) C^ )T is convex and max-convex, since 

C *,r- = PK X G int(M") | SijXj < rxi}, (11) 

where the sets whose intersection is taken are convex and max-convex. P^, r inherits these 
properties coordinate section of r . 

(ii) follows from Proposition 12.21 applied to S; (iii) follows from pj] Theorem 1 applied 
to S. □ 

We remark that sets that are both convex and max-convex have appeared under 
various names like Kleene cones [28], polytropes [IS] , or zones [23] . 
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3 Globally Optimal Solutions 



The application of the max algebra to the AHP is motivated in [mi E] by the follow- 
ing considerations. First, it is observed that, for an 57?- matrix A, vectors in the set Ca 
minimise the function ([1]) and hence the relative error. Based on this observation, these 
vectors are used to construct transitive matrices to obtain an overall ranking of the alter- 
natives in the decision process. In light of the properties of C A , this is justified by the fact 
that the transitive matrices constructed in this way are closest to the original 57?-matrix 
A in the sense of the relative error. 

Thus, the approach to construct a ranking vector for a single S'.ff-matrix taken in 
[TOl [TT] amounts to solving the following optimisation problem. 

min \e A (x)}. (12) 

xGint(R^) 

In this and the following section, we are concerned with extending the above approach to 
the general AHP with n alternatives and m criteria. 

Formally, we are given m 572-matrices; one for each criterion. Let \I/ in ([6]) denote 
the set of these matrices. For each Aj G there is an error function e^. : int(R") — > 
R+ defined as in (CQ). In contrast to the approach taken in the classical AHP, we view 
the construction of a ranking vector for the m criteria as a multi-objective optimisation 
problem for the error functions e^ 4 , 1 < i < m. 

To begin with, we seek a vector that simultaneously minimizes all of the functions e Ai • 
Such a vector is said to be a globally optimal solution for the multi-objective optimisation 
problem. 

Note that for each Ai G \&, the set of vectors that minimise e Ai : int(R") — > R + is 
precisely [H] : formally, 

C Ai = {x G int(R") | e Ai (x) = min e Ai (w)}, z = l,2,...,m. (13) 

Hence, the problem of finding a vector x G int(R") that simultaneously minimises all the 
error functions amounts to determining when 

m 

i=i 

Equivalently, x simultaneously minimises all the error functions if and only if it is a 
common subeigenvector of for all i G {1,2, ...,m}. The remainder of this section is 
divided into two parts: we first consider the existence of common subeigenvectors for 
arbitrary nonnegative matrices in the next subsection; we then specialise to sets of SR- 
matrices and globally optimal solutions. 
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3.1 Common Max-algebraic Subeigenvectors of Nonnegative Ma- 
trices 

First of all, we consider the general problem of finding a common subeigenvector for a set 
of nonnegative matrices (not necessarily 572-matrices) . Our results are clearly related to 
the work in [19] concerning the intersection of eigencones of commuting matrices over the 
max and nonnegative algebra. 

In the next result, we adopt the notation Af = ^j-y for 1 < i < m and, in an abuse 

m 

of notation, S = ® Aj. 

i=l 

Theorem 3.1. Consider the set ^ C M" xn in (0|). The following assertions are equivalent. 

(i) MS) = i; 

(ii) There exists some x G int(IR™) with A ® x < fi(Ai)x for all Ai G ^; 

(Hi) fi(Aj 1 <g) • • • (g) Aj p ) < ^(AjA ■ ■ ■ fi(Aj p ) where 1 < jk < m for 1 < k < p. (We say 
that n is submultiplicative on 

Proof. (i)=^(ii): First, assume that fi(S) = 1. Then, there exists x G int(R?) such that 
S <S> x < x. Thus, A\®x<S®x<x for all Ai G Hence, Ai <8> x < fi(Ai)x for all i. 

(ii) =^(iii): Suppose that there exists some x G int(R+) with Ai <8> x < /i(Ai)x for all 
Ai G \P. Pick some ip G ^ such that ip — Aj x £§> Aj 2 <g> • • • (g) Aj where 1 < jk < m for 1 < 
k < p. Then, 

ip ® x = A h ® A j2 <g> • • • ® A ip g> x 

< //(A^JA^ <g> A i2 <g> • • • g> A iv _ x ® x 

< ^{A jp )^{A jp _A ■ ■■fi(A j2 )A jl <g> x 

Writing r = ^(A^) fj,(Aj 2 ) ■ ■ ■ fi(Aj ), we see that x G C^ r from the definition (j2J). Hence, 
Cip,r 0- Point (i) in Proposition 12.11 implies r > fJ,(ip). Thus, fj,(AjAfj,(Aj 2 ) ■ ■ • fi(Aj ) > 
//(A?! <g) A,- a <g) • • • <g ). Note that we essentially used (@J. 

(iii) =^(i): Consider the set of normalised matrices 

if = {ii,i 2 ,...,i m } 
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where fi(Ai) = 1 for all i £ {1,2, ...,m}. Pick some ip £ As fi is submultiplicative 
on it is also submultiplicative on Thus, we have < 1. As this is true for any 
■0 £ it follows that 

max /i(?/>) < 1. 

Taking the p th root and lim sup of both sides, we see that 

Theorem 12.11 (ii) then implies that /i(S') < 1. Furthermore, since fi(S) > /i(Aj) = 1 for 
all i £ {1,2,..., m} we obtain fi(S) = 1. □ 

Note that the equivalence (i)-v4>(ii) can be regarded as a special case of Hershkowitz- 
Schneider |T7] Theorem 2.5, see also [6] Theorem 3.5 for an extension. In these works, 
the problem of simultaneous nonnegative matrix scaling is considered; this amounts to 
finding a diagonal matrix X such that XAkX^ 1 < for k — 1, . . . , m. For our case, 
take Bk = l n xn (the all-ones matrix) and impose n(Ak) = 1. However, condition (iii) 
does not appear in [17] or [6]. 

In the terminology of Butkovic et al. [51 [28], there exists a simultaneous visuali- 
sation of all of the matrices in meaning that X~ x AiX < fi(Ai) ■ l nxn for all i, and 
in particular (X" 1 AiX)ki = fi(Ai) for all % and (fc, V) £ E c (Ai). The following result 
for general nonnegative matrices will be useful in the next subsection to clarify the rela- 
tionship between commutativity and the existence of globally optimal solutions for 3x3 
STi-matrices. 

Proposition 3.1. Let A, B £ M+ Xn have /i(A) = /2(B) = 1. Iffi(A®B) = 1, then 

(i) aij = bij for all edges £ E C (A) n E C (B); 

(ii) aijbji = 1 for £ E°(A) and £ E C (B). 

Proof. As, n(A) = 11(B) = 1, it follows that A = A and B = B. Thus, S = A © B. From 
the assumption, we obtain 

n(s) = 1. 

It now follows from Theorem 13. II that there exists some x > with A®x < x, B®x<x. 
Let X = diag(x) and consider the diagonally scaled matrices 

X^AX, X~ l BX. 
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From the choice of X it is immediate that 



rtX<l nxni X- X M<W (14) 

Furthermore, from well-known facts about diagonal scaling (see Proposition 2.10 of [28J), 
it follows that 

1 = n(A) = niX^AX), 1 = fi(B) = n(X- l BX) (15) 

and that 

E C (A) = E C (X~ 1 AX), E C (B) = E C (X~ 1 BX). (16) 

We prove (i): Let G E C (A) H E C (B) be given. It follows from flT§D that 

is also a critical edge in the digraphs of X~ x AX and X~ l BX. ( fT4l) and ( Tl5l) now imply 
that 

Xi Xi 

Hence 

Xj 

and this completes the proof. 

We prove (ii): Let G E C (A) and (j,z) G E C (B). It follows from ([16]) that (i,j) 

is also a critical edge in the digraph of X~ x AX and (j, i) is a critical edge in the digraph 
oiX~ x BX. Then 

&ijXj bjjX{ . 

-Li 



and hence 



□ 



We next recall the following result, which was established in [19] and shows that 
commutativity is a sufficient condition for the existence of a common eigenvector for 
irreducible matrices. 

Proposition 3.2 ([IS]). Consider the set * C M+ Xri in (TJ|). ^sswme that each A { G * is 
irreducible and moreover that 

Ai <g> Aj = Aj <g> Ai for 1 <i,j <m. (17) 

Then there exists some x G int(M!L) with Ai® x = ^{AAx for 1 < i < m. 
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The next corollary is an immediate consequence of Proposition 13.21 and the fact, which 
we recalled in Section [TJ that for an irreducible matrix A, the set Ca is the subeigencone 
V*(A), which coincides with the eigencone V(A*). 

Corollary 3.1. Consider the set \l> C M" xn in Assume that each Ai G ^ is irreducible 
and moreover that 

A* g> A* = A* <g> A* for \<i,j<m. (18) 
Then there exists some x G int(IR™ ) with A^® x < fi(Ai)x for 1 < i < m. 

Note that flU]) implies (jlgjl. 



3.2 ^-matrices and Globally Optimal Solutions 

In the remainder of this section, we will only focus on ^-matrices. We first present the 



following corollary of Theorem 13.11 which develops the concept of simultaneous visualiza- 
tion for S'iJ-matrices. Before stating the corollary, define the anticritical graph of an 
SR- matrix to consist of the edges E C (A) given by: 



z,j)eE c (A)^(j, i )eE^\A) 



C i 



(19) 



Corollary 3.2. Consider the set \P C M^ xn in Assume that each A, G $ is an 
SR-matrix. If any of the equivalent statements of Theorem \3.1\ holds, then there exists 
some x G m£(R") such that for X = diag(x) we have 



fT^Ai) ■ lnxn < X^AX < n(Ai) ■ l r 



(20) 



In particular, 



(k,l) G E c (A i ) & (X-'A^m = ^ 
(kj) G E°{Ai) & (X- l AiX) kl = ^-\Ai 



(21) 



Proof. The right-hand side inequality of (1201) is the same as Theorem 13.11 (ii). For the 
remaining left-hand side inequality of (120]) we observe that x~ x a{jXj < n(A) is equivalent 
to xj a^~: Xi > n~ l (A). Then we apply a~- 



We next show that two distinct 572-matrices A, B in 



d2x2 



□ 



cannot have a common 



subeigenvector. Let 



A 



1 a 




1 b ' 


,B = 




l/a 1 




1/b 1 
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and assume that A ^ B. Clearly, fi(A) = fi{B) = 1 and S = A © B. If a > b, then 
1/a < 1/b and fi(S) = a/b > 1. If b > a, then 1/b < 1/a and fi(S) = b/a > 1. In both 
cases, fi(S) ^ 1. Hence by Theorem 13.11 A and 5 do not have a common subeigenvector. 

Proposition 13.21 shows that commuting irreducible matrices possess a common max 
eigenvector. We now show that for 3x3 ^-matrices, commutativity is both necessary 
and sufficient for the existence of a common subeigenvector. 

Remark 3.1. For an SR-matrix A e M^ x3 ; it is immediate that all cycle products of 
length one and two in D(A) are equal to 1. Further, there are two possible cycle products 
of length 3 in D(A): 0,12023031 and 013032021- As A is an SR-matrix, it follows that 

1 

012023031 — 

Ol3 a 32 a 21 

and hence at least one of the above products must be greater than or equal to 1. Since 
fi(A) > 1, one of the cycles of length three is critical, and the other cycle is anticritical. 
Thus, N C (A) = N(A) and A c is irreducible. Hence, it follows from Proposition ^. %\ that 
A has a unique subeigenvector up to a scalar multiple in Ca which is its max eigenvector. 
Observe that each edge (i,j) with i 7^ j belongs either to the critical or to the anticritical 
graph. 

Our next result characterises when two ^-matrices in M^ x3 have a common subeigen- 
vector. 

Theorem 3.2. Let a set {A, B} C M 3 ^ 3 of SR-matrices be given. Write S = A@B. The 
following are equivalent. 

(i) /i(5) = l; 

(ii) A and B commute; 

(Hi) There exists a vector x 6 int(W£) with A <g> x < fi(A)x, B <g> x < fi(B)x. 

Proof. The equivalence of (i) and (iii) follows immediately from Theorem 13.11 so we will 
show that (i) and (ii) are also equivalent. 

(ii) =>- (i) follows immediately from Proposition 13.21 We prove (iii)=^(ii). First note 
that it follows from Remark 13.11 that for distinct i,j,k, the edges (j,k) are either 

both critical or both anti-critical for A. The same is true of B. Calculating X _1 AX and 
X~ 1 BX where X = diag(x), it follows from Theorem 13. 1[ Corollary 13.21 and the identities 
fji{A)fji{B) = fi(B)fjL(A), fi(A)^ 1 {B) = n~ l {B)n{A) that 

dijbjk = bijQjk (22) 
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for any distinct i,j,k. It now follows from f )22|) that for i ^ j 

(A © B)ij = Clubij © a ijbjj ffi dikbkj 

= 6jj © Ojj © bikdkj 
= (B® A)ij 

where k ^ i, k ^ j . Rewriting (I2"2"j) as OjAfc = b^a^ it follows readily that fr^a^ = Oyfyj 
and ajfcfefci = h^a^. It now follows that for 1 < % < 3, 

(A © = Ojjfejj © ciijbji © cii^bki 
= buciii © bijCiji © bikCiki 
= (B® A) u 

Thus, A© J B = J B©A as claimed. □ 

It is now straightforward to extend the above result to an arbitrary finite set of SR- 
matrices in M+ x3 . 

Theorem 3.3. Let a set {Ai,...,A m } C M^ x3 of SR-matrices be given. Write S = 
Ai © • • • © A m . T/ie following are equivalent. 

(i) fx(S) = 1; 

(mJ Aj © = Aj © Aj /or a// z, j; 

(Hi) There exists a vector x G m£(R") ratt Aj © x < /i(Aj)x for all i. 

Proof. As above, the equivalence of (i) and (hi) follows immediately from Theorem 13.11 
and (ii) =>- (i) follows immediately from Proposition 13.21 To show that (i) =>- (ii), suppose 
fi(S) = 1. Then it follows that for all i,j, 

Ai © Aj S 

and hence that /x(Aj © Aj) < 1. As /i(Aj © A,-) > 1, it is immediate that 

/j(Ai ®Aj) = l 

for all i,j in {1, . . . , m}. It follows immediately from Theorem 13.21 that 

Ai © Aj = Aj © A { 

for 1 < i,j < m as claimed. □ 
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We note with the following example that commutativity is not a necessary condition 
for 4x4 572-matrices to possess a common subeigenvector. 

Example 3.1. Consider the SR-matrices given by 





1 


8 


1/4 


7 




1 


4 


5 


9 


A = 


1/8 


1 


6 


1/4 


B = 


1/4 


1 


1/8 


9 




4 


1/6 


1 


4 




1/5 


8 


1 


V 




.1/7 


4 


1/4 


1 




. V 9 


1/9 


8 


1 



where n(S) = 1. #ere ; a; = 1 0.721 0.693 0.667 
However, it can be readily verified that A® B ^ B £g> A. 



is a common subeigenvector. 



4 Min-max Optimal Points and the Generalised Spec- 
tral Radius 

In general, it will not be possible to find a single vector x that is globally optimal for the 
set \& of ^-matrices given by ([6]). With this in mind, in this short section we consider a 
different notion of optimal solution for the multiple objective functions : int(M™) — > 
1 < i < m. In fact, we consider the following optimisation problem. 



min max CaAx) . (23) 

xeint(R™) \l<i<m 'J 

In words, we are seeking a weight vector that minimises the maximal relative error where 
the maximum is taken over the m criteria (57?- matrices). 

Corollary 12.11 has the following interpretation in terms of the optimisation problem 
given in (|23|) . 

Proposition 4.1. Consider the set \1/ given by (0|). Then: 



(%) u(fy) = min ( max .(x x 

xgint(R") \l<i<m lV 

(ii) x solves (G2|) if and only if x G Cy. 
Proof. Corollary 12 . 1 1 shows that there exists some x G int(lR™) with 

max eAi(x) < r 

l<i<m 

if and only if r > fi{^)- (i) follows from this observation. The result of (ii) is then 
immediate from the definition of C^. □ 
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5 Pareto Optimality and the AHP 



Thus far, we have considered two different approaches to the multi-objective optimisation 
problem associated with the AHP. In this section we turn our attention to what is arguably 
the most common framework adopted in multi-objective optimisation: Pareto Optimality 
[U [22]. As above, we are concerned with the existence of optimal points for the set of 
objective functions e Ai , for 1 < i < m associated with the set \& ([6]) of SR- matrices. We 
first recall the notion of weak Pareto optimality. 

Definition 5.1 ([22]). w £ int(M") is said to be a weak Pareto optimal point for the 
functions e Ai '■ int(WV) — > < i < m) if there does not exist x £ int(R") such that 

e Ai (x) < e Ai {w) 

for all i = 1,2, m. 

The next lemma shows that every point in the set is a weak Pareto optimal point 
for e Al , . . .,e Am . 

Lemma 5.1. Let ^ C M™ xn be given by (T2|). Any w £ is a weak Pareto optimal point 
for e Al , . . .,e Am . 

Proof. Let w £ Cq, be given. Then e Ai (w) < fi(^) for 1 < i < m. If there exists some 
x £ int(IR™) such that e Ai (x) < e Ai (w) for 1 < i < m, then for this x 

e Ai (x) 

for 1 < i < m. This contradicts Proposition 14.11 □ 

We next recall the usual definition of a Pareto optimal point. 

Definition 5.2 ([22]). w £ int(lR™) is said to be a Pareto optimal point for the functions 
e Ai : int(R") — > K + (l <i<m) if e Ai (x) < e Ai {w) for 1 < i < m implies e Ai (x) = e Ai {w) 
for all 1 < i < m. 

We later show that the multi-objective optimisation problem associated with the AHP 
always admits a Pareto optimal point. We first present some simple facts concerning such 
points. 

Theorem 5.1. Let \I> C M" xn be given by (TJj. Then: 
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(i) If w G is unique up to a scalar multiple, then it is a Pareto optimal point for 

&Ai j • • • > e A m ; 

(ii) If w G Ca; is unique up to a scalar multiple for some % G {1,2, ...,m}, then it is a 
Pareto optimal point for ca x , • • • , ^A m ■ 

Proof. Observe that both conditions imply /t(\l/) > 0. 

(i) Assume that w G is unique up to a scalar multiple. Pick some x G int(M") such 
that eAi( x ) — e Ai{w) for all i. Then, e^.(a;) < A(^0 f° r a ^ * which implies that x G C*. 
Thus, x = aw for some a G R+. Hence, e^x) = e^io) for all i and w is a Pareto 
optimal point. 

(ii) Assume that for some 6 w 6 Ca; is unique up to a scalar multiple. Suppose 
x G int(IR™) is such that e^. (V) < eWw) for all 1 < j < m. In particular, a; G Ca,, and 
this implies that x = aw for some a G M. Further, it is immediate that for any other 
Aj G ^ {i 7^ j), we have e^fx) = ^ (w). Thus, w is a Pareto optimal point. □ 

By Proposition 12.21 condition (i) is equivalent to N C (S) = N(S) and S to be irre- 
ducible, and condition (ii) is equivalent to N c (Ai) = N(Ai) and Af to be irreducible for 
some i. 

Corollary 5.1. Let the set \l/ C M™ xn given by ftfity consist of SR-matrices. Forn G {2, 3}, 
any w G Cj^il < i < m) is a Pareto optimal point for , • • • , &A m ■ 

Proof. Notice that from Remark 13.11 for 3x3 case, there exists a unique subeigenvector 
up to a scalar multiple in each for 1 < % < m. This is also true for the 2x2 case 
because N C (A) = N(A) and A° is irreducible. The result directly follows from (ii) in 
Theorem O □ 

The following example demonstrates point (i) in Theorem 15.11 
Example 5.1. Consider the following matrices given by 



A 



1 9 1/4 2 

1/9 16 3 

4 1/6 1 1/4 

1/2 1/3 4 1 



B 



1 1/2 4 1/8 

2 13 2 
1/4 1/3 1 5 

8 1/2 1/5 1 



S matrix is obtained as follows 



1 



9 



4 2 



S 



2 



1 



6 3 



4 1/3 15 
8 1/2 4 1 
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where N (S) = N(S) and S is irreducible. From Proposition we have a unique 

Figure U\ below 



1 T 



1 0.758 0.861 1.174 



vector (up to a scalar multiple) in Cq,: w - 

represents the values of e&{x) and cb{x) at w and some points in Ca, Cb and mt(L + , 
Remark that Pareto optimality is observed at w G where ca{w) = esiw) = 6.817. 
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Figure 1: Refer to Example 15.11 



Our objective in the remainder of this section is to show that the multi-objective 
optimisation problem associated with the AHP always admits a Pareto optimal solution. 
We first recall the following general result giving a sufficient condition for the existence 
of a Pareto optimal point with respect to a set £ C M". Essentially, this is a direct 
application of the fact that a continuous function on a compact set always attains its 
minimum. 

Theorem 5.2 (|4J). Let E C MJ 1 be nonempty and compact. Let a set of continuous 
functions {fi, . . . , f m } be given where 

f r .E^R + 

for 1 < % < m. There exists w G E such that x G E, fi(x) < fi(w) for 1 < i < m implies 
fi(x) = fi{w) for 1 <i <m. 
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This result follows from elementary real analysis and the observation that if w min- 

m 

imises the (continuous) weighted sum ^ &ifi{x) where a, > for 1 < i < m, then w 

i=i 

must be Pareto optimal for the functions fi, ■ ■ ■ , f m . 

A point w satisfying the conclusion of the Theorem 15.21 is said to be Pareto optimal 
for {/i, . . . , f m } with respect to E. Thus, for any multi-objective optimisation problem 
with continuous objective functions defined on a compact set, there exists a point that is 
Pareto optimal with respect to the given set. 

To apply the above result to the AHP, we first note that for a set \& of 572-matrices, T>^ 
is compact by Theorem 12.31 (recall that STi-matrices are positive and hence irreducible). 
We now show that any point in that is Pareto optimal with respect to T>^ is also 
Pareto optimal with respect to the set int(R"). 

Although we don't specifically use the SR property in the following, we assume \I/ to 
consist of 57?-matrices as we are primarily interested in the AHP application. Instead, 
we could just assume that Ai are irreducible implying that T>^ is compact and a Pareto 
optimal point exists. 

Lemma 5.2. Consider the set \l/ in (GJ) and assume that Ai is an SR-matrix for 1 < i < m. 
Let w be a Pareto optimal point for {e^, • • • , eA m } with respect to V^,. Then, w is also 
Pareto optimal for e^, . . . , CA m with respect to int{W\_). 

Proof. Assume that w G T)^ is a Pareto optimal point with respect to V^. Suppose 
x G int(lR™)\C*. Then from the definition of (jHJ), it follows that 

&Ai (x) > /H 1 ^) f° r some io- (24) 

As w G caXw) < fi{^f) for 1 < i < m. It follows immediately from (12411 that for any 
x ^ it cannot happen that e^ipc) < eAX w ) for 1 < i < m. 
Let x in Cq, be such that 

e Ai(x) < eAi(w) for 1 < i < m. 

As eAi(Xx) = eAi(x) for all A > 0, 1 < i < m, and w is Pareto optimal with respect to 
Uy, it follows that e^x) = e^(w) for 1 < % < m. □ 

Our next step is to show that there exists a point x G T>* that is Pareto optimal with 
respect to T>^. 

Proposition 5.1. Consider the set ^ in |6]) and assume that Ai is an SR-matrix for 
1 < % < m. There exists x G T>^ that is Pareto optimal for e^, . . . , CA m with respect to 
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Proof. First note that T>^, ^ since /t(^) > 0. Theorem 12.31 shows that V>^ is compact. 
Furthermore, for any irreducible matrix A, the function : — > R + is a composition 
of continuous functions and hence continuous on a compact set. Theorem 15.21 implies that 
there exists w in T>^ that is Pareto optimal with respect to Vq,. □ 

Combining Proposition 15.11 with Lemma I5.2[ we immediately obtain the following 
result. 

Corollary 5.2. Consider the set \I> in (G|) and assume that Ai is an SR-matrix for 1 < i < 
m. There exists x G Dq, that is Pareto optimal for e Al , . . . , e^ m with respect to int(W±). 

Corollary 15.21 means that there exists a vector x of positive weights that is simulta- 
neously Pareto optimal and also optimal in the min-max sense of Section H] for the error 
functions e^, 1 < i < m. 

Finally, to illustrate the above results, we revisit Example 11.11 

Example 5.2 (Example 11.11 Revisited). Let C, Ax, . . . , A 5 be as in Example \l.l\ Taking 
on, 1 < i < 5 to be the ith entry of the max eigenvector of C , normalised so that ai = 1, 
we apply Theorem \5.S\ to compute Pareto optimal solutions in the set T>^ by minimising 
the weighted sum 

m 

^2a>ie A .(x) 

i=l 

using the MATLAB function fminsearch. 

Observe that ji{S) = 4.985, so there is no common subeigenvector in this case. Next, 
we calculate the max eigenvector of C : 

1 

1.495 
2.236 
3.344 
0.897 

We find that there are multiple Pareto optimal points giving at least two possible distinct 
rankings: 1 > 3 > 4 > 2 and 1 > 3 > 2 > 4. Notice that first ranking scheme is the same 
as the one obtained from the classical method used in Example The second ranking 
scheme is also reasonable, since if we analyse the local rankings associated with the set of 
SR-matrices in detail, we see that 2 > 4 for A\, A 2 and A 5 . In particular, 2 is preferred 
to all other alternatives for A 2 ■ 
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6 Conclusions and Future Work 



Building on the work of Eisner and van den Driessche PJJJ ELT] , we have considered a max- 
algebraic approach to the multi-criteria AHP within the framework of multi-objective 
optimisation. Papers [TOT ITT] characterise the max eigenvectors and subeigenvectors of 
a single SR-maXhx as solving an optimisation problem with a single objective. We have 
extended this work to the multi-criteria AHP by directly considering several natural ex- 
tensions of this basic optimisation problem to the multiple objective case. Specifically, we 
have presented results concerning the existence of: globally optimal solutions; min-max 
optimal solutions; Pareto optimal solutions. The principal contribution of the paper is 
to draw attention to this max-algebraic perspective on the multi-criteria AHP, with the 
main results in this direction being: establishing the connection between the generalised 
spectral radius and min-max optimal solutions (Proposition 14. IP ; proving the existence of 
Pareto optimal solutions and showing that it is possible to simultaneously solve the Pareto 
and min-max optimisation problems (Proposition 15.11 and Corollary 15 . 2 j) . We have also 
related the existence of globally optimal solutions to the existence of common subeigen- 
vectors and highlighted connections between this question and commutativity (Theorem 
HO}. 

Buket Benek Gursoy and Oliver Mason acknowledge the support of Irish Higher Educa- 
tional Authority (HEA) PRTLI Network Mathematics Grant; Sergei Sergeev is supported 
by EPSRC Grant RRAH15735, RFBR-CNRS Grant 11-01-93106 and RFBR Grant 12- 
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